High cooperativity coupling of electron-spin ensembles to superconducting cavities 
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Electron spins in solids are promising candidates for quantum memories for superconducting 
qubits because they can have long coherence times, large collective couplings, and many quantum 
bits can be encoded into the spin-waves of a single ensemble. We demonstrate the coupling of elec- 
tron spin ensembles to a superconducting transmission-line resonator at coupling strengths greatly 
exceeding the cavity decay rate and comparable to spin linewidth. We also use the enhanced cou- 
pling afforded by the small cross-section of the transmission line to perform broadband spectroscopy 
of ruby at millikelvin temperatures at low powers. In addition, we observe hyperfine structure in 
diamond PI centers and time domain saturation-relaxation of the spins. 



An eventual quantum computer, like its classical ana- 
log, will make use of a variety of physical systems special- 
ized for different tasks. Just as a classical computer uses 
charge-based transistors for fast processing and magnetic 
hard drives for long term information storage, a quantum 
computer might use superconducting qubits for process- 
ing [1 and ensembles of electron spins as quantum memo- 
ries [2, 3 , linked by single microwave photons. Although 
other microscopic systems have been proposed for use in 
a hybrid architecture [IHZ], electron spins complement 
superconducting qubits particularly well. They feature 
similar transition frequencies, do not require trapping, 
and can be packed densely. Furthermore, a single ensem- 
ble could be used to store many qubits using holographic 
encoding techniques [3 demonstrated classically for nu- 
clear [8 and electron [9 spins. 

In this Letter, we demonstrate the first step toward re- 
alizing a solid-state quantum memory: coupling an elec- 
tron spin ensemble to an on-chip superconducting cavity 
at powers corresponding to a single cavity photon. We 
observe megahertz spin-photon interaction strengths in 
both ruby Cr^+ spins and N substitution (PI) centers 
in diamond. A parallel effort by Kubo, et. al.[? ] sees 
similar coupling to nitrogen vacancy (NV) centers in dia- 
mond. In doing so we develop a platform for the study of 
electron spin resonance (ESR) physics in picoliter mode 
volumes, millikelvin temperatures, and attowatt powers. 
Finally, we perform time-resolved saturation/relaxation 
measurements of the PI centers, a precursor to full pulsed 
control of the system. 

ESR studies the microwave response of electron spins 
at their resonant frequency in a magnetic field. Samples 
are conventionally placed inside a 3D high quality-factor 
(Q) cavity which enhances the sensitivity by confining 
photons with the spins and extending the interaction 
time [l^. In this work, several ID cavities are capaci- 
tively coupled to a common feedline on a sapphire chip. 
We place the spins within the mode volume by fabri- 
cating the device on doped sapphire (ruby - Fig. [l^), at- 
taching a substrate on top of an existing device (diamond 



- Fig. [T]3), or simply spin-coating the surface (DPPH - 
not shown). The single spin-photon coupling is given by 
^s/27r = mo{jj.oUj/2hVcy^'^^ where mo is the spin's mag- 
netic dipole moment. The single spin coupling can be 
enhanced by using small mode volume cavities [FH] |T2] . 
Using a superconducting cavity allows use of extreme as- 
pect ratios while maintaining high Q. In this case by 
using a ID {Vc ~ d'^X) rather than a 3D {Vc ~ A^) cavity, 
we increase the coupling by the ratio of the wavelength to 
the cavity width, X/d ^ 1000. Using an ensemble the ef- 
fective collective coupling is ^s,eff ^ M^^'^Qs- Employing 
a large number of spins M - 10^^-^^ enables megahertz 
interaction strengths. Because the resonators are made 
from superconducting thin films, they can maintain a 
high Q despite their extreme aspect ratio and ~ 200 mT 
applied in-plane magnetic fields [13 . 

Transmission through the ruby chip feedline is plotted 
in Fig. [T]^ from 10 to 14.5 GHz as a function of the in- 
plane magnetic field. The most striking features are the 
transitions of the spin-3/2 Cr^+ which are visible over a 
broad frequency range. 

They can be identified by diagonalizing the ruby 
Hamiltonian [M] 

Hr = -mo,rB-S-D{Sl-^), (1) 

with mo,r/27r = 27.811 MHz/mT, S is the spin-3/2 
operator, z is defined along the c axis, and 2D/2t: = 
11.46 GHz. In addition to the Zeeman term there is a 
large crystal field in ruby that separates the levels into 
two doublets {S^ = ±3/2 and = ±1/2). The mea- 
sured spectroscopic lines correspond to transitions be- 
tween the doublets. When a magnetic field is applied at 
angle with respect to the c axis the levels within each 
doublet are hybridized giving the observed curvature in 
the transitions. Because of the relatively high doping 
~ ]^q19-20 (3j^-3 because the spins efficiently fill the 
mode volume of the feedline, absorption at the transition 
frequencies of the Cr^+ spins is visible even without the 
aid of a cavity. From the 5% absorption dips we estimate 
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FIG. 1. Cavity-spin devices and broadband spectroscopy of 
ruby, a) Five A/4 resonators with fundamental frequencies 
from 4 to 8 GHz coupled to a common feedline fabricated on 
a ruby substrate. The coupling Q is determined by the length 
of the elbow coupler running parallel to the feedline. b) A 
piece of synthetic diamond glued on top of a similar device on 
an undoped c-plane sapphire substrate, c) The transmission 
spectrum of the ruby device when a magnetic field is applied 
parallel to the Nb thin film. The dashed curves indicate ruby 
transitions. Fit uses only two parameters, the magnet cur- 
rent to field ratio and the angle between the magnetic field 
and crystal c axis which is found to be 64.6°. To compensate 
for frequency dependent attenuation in the microwave lines, 
each point is divided by the mean of all points with the same 
frequency (row). The three broad horizontal lines are reso- 
nances with the copper sample holder, whereas the narrow 
line at ^ 14.35 GHz is the second mode of the longest super- 
conducting transmission line resonator. A higher resolution 
scan of the area within the red box is shown in Fig. |2^. 



that there are approximately 10^^ spins interacting with 
the feedline. 

In addition to the spectrum measured via the feed- 
line, several magnetic field insensitive resonances are ev- 
ident in Fig. [T]^. The bottom three broad modes are 
resonances in the sample holder, while the narrow mode 
at 14.35 GHz is the superconducting transmission line 
cavity. Both the modes of the copper sample holder and 
the superconducting cavities can be used to manipulate 
spins, but we focus on the on-chip cavities which are 
higher Q and have typically been used in qubit exper- 
iments [l]. 
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FIG. 2. a) Transmission spectrum as a function of frequency 
and magnetic field near resonance with the second mode of 
one of the superconducting cavities. The nearly vertical fea- 
ture in the transmission spectrum at 60 mT is due to the 
1^4 transition in the spins coupled to the feed line. An 
avoided crossing with the superconducting cavity is visible 
indicating coupling to M = 4 x 10^^ spins. The abrupt jumps 
in the cavity frequency are thought to be due to flux penetrat- 
ing the superconducting film causing changes in the effective 
inductance of the cavity, b) Fit using Eq. [5] to the Q of the 
cavity, which is partially damped by the spins. 

In a higher resolution scan (Fig. [2|i) an avoided cross- 
ing is present where the 1^4 transition approaches 
the cavity at 14.35 GHz. Also visible are the interact- 
ing spins within the feedline which are unaffected by the 
cavity-spin coupling. The cavity line in Fig. [2^ is repelled 
by more than its linewidth and on resonance it is damped 
primarily by the spins. Because of flux jumps (due to a 
small component of the field perpendicular to the super- 
conducting film), it is difficult to precisely fit frequency 
versus field. However, the Q of the cavity is relatively 
unaffected by these jumps, and can be used to extract 
the same information. 

The cavity Q is given by 

^ 2 I 2 

where both the spins and the cavity are modeled as 
single-mode harmonic oscillators with detuning A, cav- 
ity frequency Co'^/27r = 14.35 GHz, and spin frequency 
degenerate at B=60 mT with the spin resonance tuning 
at rate meff/27r = 52.4 MHz/mT. The cavity linewidth 
K/27r = uJr/27rQ = 1.3 MHz is independently measured 
away from resonance. The collective coupling strength 
and spin decoherence rate are determined from the fit to 
be ^s,eff/27r = 38 MHz and 72/27r = 96 MHz. The latter 
is probably primarily due to broadening by strong hyper- 
fine interactions with ^^Al nuclear spins The cou- 
pling is much larger than the decoherence rates measured 
in both the cavity and typical superconducting qubits [l] . 
A dimensionless measure of the coupling strength in cav- 
ity QED [? ] is the cooperativity C = g'^ / ^ 11-5. 
Since C > 1 the coupling is strong in the sense that 
at resonance nearly every photon entering the cavity is 
coherently transferred into the spins. In order to make 
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FIG. 3. Experimental (a) and analytically calculated (b) 
transmission spectrum of cavity- diamond spin system. In 
both plots, 5 spin resonances due to the anisotropic hyper- 
fine splitting are observable (see text for details) . The dashed 
lines in (b) refer to bias fields used in Fig. [4] 
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retrieval possible a spin system with smaller linewidth is 
required. 

Using the diamond sample we are able to observe 
hyperfine splittings and perform time-domain satura- 
tion/relaxation experiments. ESR active PI centers [16] 
are visible via absorption spectroscopy (Fig. [3|, but we 
did not observe resonance with NV centers, which have 
been studied in a parallel publication [? ]. 

The observed splittings can be understood by consid- 
ering the transitions of the Hamiltonian [16] 

= -mo^dB - S ^ AS ■ I, (3) 

where mo,d/27r = 28.04 MHz/mT and the hyperfine cou- 
pling tensor A/27r = (81.33,81.33,114.03) MHz. Here, 
the z direction corresponds to the diamond 111 axis, S 
are the electron spin- 1/2 operators, and / are the nuclear 
spin-1 operators. This describes a nitrogen atom with an 
extra, nearly-free electron, and a hyperfine spectrum due 
to the 1 = 1 nuclear spin of ^^N. This splits the line into 
three: m/ = —1,0,1. The outer two lines are further 
split by anisotropy in the hyperfine coupling depending 
on which carbon the nitrogen substitutional's electron 
prefers. Selection rules suppress transitions which change 
m/, giving a total of five lines shown in Fig. |3]3, and ob- 
served in Fig. |3^. 

The model in Fig. [sJd treats each of the twelve possible 
transitions (three hyperfine transitions and anisotropics 
due to the four bond angles) as an independent harmonic 
oscillator with each oscillator uniformly coupled to the 
cavity. When aligned along one of the bond angles de- 
generacies reduce the spectrum to the five lines visible 
in Fig. |3| The effective coupling strengths used were 5 
MHz for the central peak, 3.5 MHz for the satellite peaks, 
along with 30 MHz linewidths for each of the transitions. 
The large linewidth is most likely due to additional para- 
magnetic impurities [17] or magnetic field inhomogeneity. 

Although this device was not designed to perform pre- 
cise rotations on the spins, it is possible to saturate the 



FIG. 4. Saturation/relaxation measurement of diamond, a) 
A saturating pulse (P ^ 100 pW) is applied to the cavity, de- 
polarizing the spins. The response of the resonator is probed 
continuously with a weak probe Prf ~ 100 fW which does not 
significantly excite the spins, b) Measurement of resonator Q 
as a function of time immediately after a 1 second saturation 
pulse. The two curves are measured at the locations of the 
magnetic field indicated in Fig. ^jp by the dashed lines where 
the spins are on (red) and off (blue) resonance with the res- 
onator, c) Measurements were taken at each frequency and 
the reconstructed cavity resonance fit at each time slice to de- 
termine the Q as a function of time. Each curve corresponds 
to a different time indicated by an arrow in (b). The off reso- 
nance (blue) curves show almost no time response and lie on 
top of each other. 



spins incoherently and observe their recovery in the time- 
domain. The data in Fig. |3^ are measured at 100 fW 
^ 100 photons in the cavity. Whereas at low powers the 
cavity is strongly damped by the spins, saturation be- 
gins at Pg = hujjj2f^/'^9s ^ 10 pW ~ 10^ photons and 
by Pg ~ 100 pW the cavity Q is restored to nearly its 
unloaded value. In Fig. [4| two experiments are shown 
at magnetic fields corresponding to the dashed lines in 
Fig. [sJd, where the spins are on (red) and off (blue) reso- 
nance with the cavity. This is done at several frequency 
points near the cavity resonance, allowing us to observe 
the evolution of the cavity lineshape with time after the 
saturation pulse. In Fig. ^ one can see that immediately 
after the saturation pulse the resonator has a high Q. As 
the spins relax the resonator becomes strongly damped 
once more on a timescale of seconds. A lifetime of ~ 1 s 
is consistent with an estimate based on the saturation 
power and the linewidth measured in Fig. [3] [18 . How- 
ever, it is difficult to extract the relaxation time quanti- 
tatively because the curve is not exponential. This could 
be attributed to the presence of other ESR centers [19], 
heating of the host lattice via the spins, or some other 
type of saturable impurity which is suspected to play a 
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role in the decay of superconducting cavities [20] . There 
appears to be a smah time dependence in the off reso- 
nance trace in Fig.[4]3. This is also consistent with broad- 
band saturable impurities ^] and might point to a 
technique by which such impurities could be temporarily 
neutralized. 

Figures [l] - [4] show that these devices can be used as 
spectrometers of electron spin resonance. The sensitivity 
of an ideal ESR spectrometer is a complicated function of 
both instrument and sample parameters, including ESR 
linewidths, spin lifetimes (which are themselves temper- 
ature and magnetic field dependent). By substituting in 
for the mode volume and saturation power the typical 
sensitivity formula in Weil [TO can be re-expressed in 
the terms of cavity QED, highlighting the fundamental 
processes in detection: 

Here, s is the total spin, r is the integration time, 
Ps = tdinh {hu / kBTs) is the polarization of the spins at 
temperature Tg, and = k^T^ /fioo is the effective num- 
ber of noise photons added by the detector. In essence 
the sensitivity is given by the number of spins required 
to give Ua scattered photons (at rate gl/n) in one in- 
tegration time. Thermal depolarization {ps)-, and inho- 
mogenous broadening (272/7) worsen sensitivity from the 
case of ideal spins at zero temperature. Conventional 
ESR instruments have Q = 10,000, ^s/27r - 0.02 Hz 
{Vc = 10 /iL), T/v = 30000 K, typical saturation pow- 
ers of Ps ~ 100 mW, and operate from 300 K down to 
1 K [To]. These typically produce an TVmin of ~ 10^ 
(polarized) spins. Our current devices have Q=5000, 
gj27r - 20 Hz (K = 10 nL), and Tn 10 K, and func- 
tion from 4 K down to 20 mK. The small mode volume 
of the ID resonators also reduces the optimal Pg allowing 
the study of samples at low temperatures where cooling 
power is limited. Under similar conditions these param- 
eters predict a spin sensitivity of A^min ~ 10^, worsened 
to about Nmin ~ 10^ due to the inhomogenous broaden- 
ing of the spins measured here. The sample contained 
approximately ~ 10^^ which were readily detected even 
without the use of standard ESR techniques such as field 
modulation. 

Through further development, future devices could 
have Q = 100, 000, = 1 kHz {Vc = 1 pL using elec- 
tron beam lithography), and = 1/2 using a quantum 
limited amplifier [21]. For these values g^/i^ ^ 100 pho- 
tons/s, and it should be possible to detect A^min ^ 1 
using a cavity and perhaps even strong coupling if su- 
perconducting qubits can be used to mediate the inter- 
actions \22j> .23^ . 

To move forward towards a quantum spin memory sev- 
eral improvements must be made. Most importantly a 
more coherent electron spin candidate is required. In 



addition, because superconducting qubits might be ad- 
versely affected by magnetic fields, an ideal candidate 
would have a ZFS such that B = Tiuo ^ k^T^ which 
could come from crystal field splittings as with the ruby 
shown here, NV~ centers, etc. Crystal field split spins 
have the advantage that they can be tuned with small in- 
plane magnetic fields allowing in-situ adjustment of the 
spin frequency and holographic encoding [3 . Other can- 
didates such as spins with large hyperfine splittings [9] 
might also be interesting as they can be insensitive to 
small magnetic field fluctuations. This comes at the ex- 
pense of not being tunable with such fields. In addition, 
they must have long coherence times, requiring control 
over line broadening due to dipolar couplings, nuclear 
spins, strain variations (for crystal field splittings) and 
inhomogeneous bias fields. New resonator designs could 
improve the static and microwave magnetic field unifor- 
mities, allowing spin-echo techniques to be applied. 

We have demonstrated coupling of large ensembles of 
electron spins to both broadband coplanar waveguide 
transmission lines and resonators. The coupling is suffi- 
ciently strong to exceed all qubit and cavity decay rates 
with large cooperativity, but is still limited by the spin 
linewidth. The system also shows promise as a general 
ESR tool allowing good sensitivity in a broadband sys- 
tem, and exquisite sensitivity in a high Q cavity. The 
small mode-volume makes it ideally suited for study- 
ing picoliter scale samples, especially 2D systems such 
as graphene or semiconducting heterostructures. Fur- 
ther applications could include maser amplification and 
single-photon microwave to optical up conversion 
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